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Few  studies  have  been  made  concerning  the  effects  of  low  cycle 
fatigue  on  the  stress/strain  field  or  the  plastic  zone  size  ahead  of 
a  crack  tip  in  a  high  temperature  environment.  Experiments  using  a 
superalloy  IN-100  at  732°C  have  shown  crack  growth  during  fatigue  is 
slower  than  the  growth  rate  during  long  periods  of  sustained  load. 

These  investigations  also  show  that  crack  growth  rates  during  hold 
periods  after  fatigue  cycling  are  dependent  on  the  cycle  frequency. 

This  study  focuses  on  fatigue  as  it  effects  the  changing  stress 
field  and  plastic  zone  ahead  of  a  crack  tip.  Finite  element  modeling 
was  accomplished  using  an  in-house  computer  program  named  VISCO,  which 
was  modified  to  incorporate  load  cycling.  A  compact  tension  specimen 
geometry  was  modeled  using  the  Bodner-Partom  viscoplastic  constitutive 
equations  to  describe  the  material  behavior.  ^  Load  spectra  consisted  of 
constant  amplitude  saw  toothed  patterns  with  non-zero  mean  load  having 
a  ratio  of  minimum  to  maximum  load  of  0.1. 

" — -)  Results  show  that  material  behavior  near  a  crack  tip  closely  cor¬ 
relates  to  uniaxial  material  response  data  under  stress  controlled 
loading.  During  cyclic  loading,  the  greatest  plastic  deformation  near 
the  crack  tip  occurs  in  the  first  load  cycle.  Computations  show  that 
IN-100  displays  significant  amounts  of  time  dependent  inelastic  behavior. 


I.  INTRODUCTION 


With  the  growing  use  of  high  performance  military  gas  turbine 
engines  and  the  escalating  cost  of  aircraft,  the  United  States  Air  Force 
places  stringent  service  limits  on  critical  airframe  and  engine  compon¬ 
ents. 

Airframe  components  are  periodically  inspected  for  flaws  and  re¬ 
turned  to  service  if  the  flaws  can  not  grow  to  critical  size  prior  to 
the  next  periodic  inspection.  Critical  engine  components,  however,  such 
as  turbine  disks,  are  removed  from  service  at  a  time  when  statistically 
1  in  1000  would  be  expected  to  have  initiated  a  crack  of  some  finite 
length  (0.03  in).  Eighty  percent  of  the  disks  have  at  least  ten  useful 
lifetimes  remaining,  though  no  attempt  is  made  to  utilize  the  statistic¬ 
ally  "failed"  disks  [2]. 

Engine  operation  occurs  in  a  complex  region  characterized  by  fre¬ 
quent  load  cycle  excursions  and  short  periods  of  sustained  loading. 

High  temperature  engine  operation  introduces  time  dependent  phenomena 
which  interact  with  varied  load  spectra  to  produce  complex  material 
behavior.  The  thrust  of  the  Air  Force  research  has  focused  on  crack 
growth  prediction  with  lesser  emphasis  on  the  interaction  of  fatigue 
cycling  and  sustained  hold  times.  If  the  effects  of  frequency  on 
fatigue  cycling  and  load  amplitude  are  fully  understood  when  coupled 
with  sustained  hold  times,  accurate  remaining  life  predictions  could  be 
made  for  components  with  subcritical  flaws.  Many  retired  parts  could  be 
kept  in  service.  Only  components  with  a  quantifiable  critical  flaw  size 
would  be  Retired-for-Cause  [14]. 


:: 


: 


Approach 

This  thesis  attempts  to  quantify  fatigue  crack  growth  and  the 
effects  of  cyclic  loading  on  IN-100,  a  superplastically  forged  Nickle- 
based  superalloy  used  in  turbine  disks  for  the  F-100  engine.  Compari¬ 
sons  between  laboratory  low-cycle  fatigue  tests  and  computer  generated 
viscoplastic  cyclic  effects  are  directed  at  the  change  of  material  par¬ 
ameters  with  frequency,  load,  and  plastic  strain.  The  coupling  of  high 
temperatures  and  stresses  causes  time  dependent  plastic  deformation  at 
discontinuities  in  the  material.  Due  to  its  independence  of  yield 
criteria  and  its  close  correlation  with  experimental  observations  of 
IN-100,  plastic  flow  near  these  discontinuities  was  modeled  using  the 
Bodner-Partom  flow  law.  In  addition,  the  Bodner-Partom  model  has  the 
capability  to  predict  behavior  produced  by  cyclic  effects,  time  depen¬ 
dent  creep,  strain  hardening,  and  plastic  deformation.  Thus,  labora¬ 
tory  low-cycle  fatigue  tests,  performed  on  compact  tension  specimens 
manufactured  from  typical  disk  forgings,  were  modeled  using  finite 
element  methods  and  the  Bodner-Partom  viscoplastic  flow  law. 

A  linear  Euler  extrapolation  scheme  was  used  to  Integrate  the 
flow  law  through  time  and  Incorporated  into  the  finite  element  model 
incorporating  the  residual  force  technique  [7]. 

VISCO,  a  computer  program  developed  by  Hinnerichs  [7]  and  used  in 
this  thesis,  uses  constant  strain  triangular  elements  and  simulates 
crack  growth  by  releasing  fixed  nodes  at  a  predetermined  time.  As  the 
node  is  released,  nodal  forces  are  redistributed  to  the  surrounding 
elements  and  the  boundary  changed  from  one  of  restraint  to  a  force 
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boundary  condition.  Costly  matrix  factorization  is  avoided  by  incor¬ 
porating  a  Gauss-Seidel  iteration  equation  solver,  in  which  pertinent 
terms  of  the  stiffness  matrix  are  the  only  ones  changed  between  time 
steps  [7].  A  cyclic  load  algorithm  was  used  to  accommodate  the  chang¬ 
ing  load  spectra.  The  fatigue  stress  spectrum  was  modeled  by  a  saw¬ 
toothed  stress-time  pattern  of  constant  amplitude  with  a  non-zero  mean 
.load  .. and.  a  load  ratio  of  .1  [14].  During  each  time  step;  stress, 
strain,  plastic  work  and  z  hardness  parameter  was  calculated  for  each 
element  in  the  model. 


II.  VISCOPLASTICITY  THEORY 


In  the  usual  sense,  plasticity  is  defined  as  permanent  deformation 
by  the  application  of  stresses  which  are  greater  than  those  stresses 
required  to  cause  yielding.  The  distortion  depends  not  only  on  the 
final  state  of  stress,  but  also  upon  the  stress  state  from  the  start  of 
yielding  [12].  Thus,  plastic  behavior  can  be  characterized  by  irre¬ 
versible  straining.  All  materials  exhibit  some  plastic  behavior,  even 
when  the  applied  loads  produce  stress  below  the  yield  stress.  Most 
plastic  deformation  below  the  yield  stress  is  small;  and  as  such,  is 
usually  neglected  [16]. 

The  total  strain  a  material  experiences  can  be  represented  by  the 


expression 


(2.1) 


£ 

where  is  the  total  strain,  e^j  is  the  elastic  (reversible)  strain 

D 

and  e-jj  is  the  plastic  (irreversible)  strain.  Viscoplastic  (non-re- 
coverable)  strains  occur  only  after  the  yield  stress  is  exceeded. 


If  one  takes  a  time  derivative  of  equation  (2.1),  a  rate  dependent 


form  for  total  strain  is  obtained 


where  ^  ^ 


(2.2) 


By  taking  a  time  derivative  of  Hooke's  law,  the  elastic  strain 
rate  c^e  can  be  related  to  the  stress  rate;  however,  it  is  necessary 
to  relate  the  plastic  strain  rate  to  stress  by  some  other  means. 

In  classical  theory,  plastic  deformation  begins  at  yield,  and  is 
dependent  on  a  yield  criterion.  General  yield  criteria  can  be  written 
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(2.3) 


f(o..)  =  K(k) 

where  f  is  some  function  of  stress,  K  is  a  material  parameter  deter¬ 
mined  experimentally  and  k  is  a  hardening  parameter.  The  yield  criteria 
as  used  herein  is  independent  of  coordinate  orientation  and  only  a 
function  of  the  three  stress  invarients  [13]. 


1 =  aij 

(2.4) 

.  1 

2  "  7  °ij°ij 

(2.5) 

.  1 

3  "  7  aijajk°ki 

(2.6) 

Through  experimental  studies  by  Bridgman,  plastic  deformation  is 
independent  of  hydrostatic  stress.  Thus,  the  yield  function  can  be 
written  as  f(02,J3)  =  K(k)  where  J2  and  J3  are  the  second  and  third 
invariants  of  the  deviatoric  stress  tensor. 


J 


„  1 


2  =  7  S1jSij 


(2.7) 


and 

Sij  *  °ij  "  I  aij°kk  (2*8 

Many  forms  of  yield  criteria  have  been  developed.  Each  predicts 
plastic  flow  under  different  conditions  and  materials.  The  Von  Mises 
yield  criterion  is  one  of  these.  It  Is  used  primarily  for  metals 
behavior  and  is  based  on  distortion  strain  energy  theory  [13].  Yield¬ 
ing  begins  when  the  distortion  strain  energy  equals  the  distortion 
strain  energy  required  to  produce  yielding  in  a  simple  uniaxial  stress 


(2.9) 


The  distortion  strain  energy  can  be  written 

ud  ■  I  i  °ijdeijdv 


or  in  terms  of  J2»  the  second  invariant  of  deviatoric  stress 


Ud  =  2G  J2 


(2.10) 


where  G  is  the  shear  modulus. 


In  explicit  terms  of  principal  stresses* the  above  expressions 


become 


J2  s  J  (°l"°2^  +  ^°2"°3^2  +  ^°ra3^ 


(2.11) 


Jd  =  T2£  (al“a2'  +  ^a2‘°3^2  +  (®r«3)‘ 


(2.12) 


For  a  uniaxial  test 


°1  *  °ys  02=03*° 


where  o  equals  the  stress  at  the  start  of  yielding. 


J2  reduces  to 


’2  "  I  °ys2 


where  oy$  is  the  uniaxial  yield  stress. 

The  distortion  strain  energy  becomes 

ii  -  1  2 

Ud  "  6G  ffys 


(2.13) 


(2.14) 


(2.15) 


Equating  the  multiaxial  distortion  strain  energy  to  that  of  the 


uniaxial  case  produces 


J2  *  K*ysC  '  (2.16) 

where  K  is  a  proportionality  constant.  Multiaxial  yielding  is  pre¬ 
dicted  when  J2  reaches  the  critical  value  at  yield  in  a  uniaxial  stress 
test. 
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Using  the  Prandtl-Reuss  relations,  which  assume  incompressibility 

p 

and  isotropy,  the  plastic  strain  e..  can  be  expressed  incrementally  as 

*  v 


(2.17) 


where  x  is  a  positive  load  history  dependent  proportionality  constant 
relating  the  material  viscosity  and  S..  are  the  components  of  the  devia- 

*  J 

p 

tone  stress  tensor.  If  given  a  form  for  dx,  the  plastic  strain  e-.  can 

*  J 

be  found  for  any  state  of  stress.  It  can  be  shown  that,  at  yield,  the 
Prandtl-Reuss  relations  imply  theVonMises  yield  criterion  [13].  The 
specific  expression  for  dx  depends  upon  the  plastic  flow  law  used. 
Bodner-Partom  Constitutive  Model 

Bodner  and  Partom  [3,4,5]  suggested  that  plastic  flow  beginsat  the 
onset  of  loading,  and  as  such,  a  continuous  flow  relationship  exists 
between  plastic  strain  and  stress.  Through  the  study  of  dislocation 
dynamics  Bodner  and  Partom  postulated  the  form  of  the  parameter  x  by 


squaring  the  Prandtl-Reuss  relation 


(2.18) 


resulting  in 


1  -P  -P  _  n  P  _  1  .Zc  c  .  .2, 

7  E1j  e1j  ■  °2  '  7  x  S1jS1J  X  J2 


(2.19) 


where  is  the  second  invariant  of  the  plastic  strain  rate,  and  J2  is 
the  second  invariant  of  the  deviatoric  stress  tensor. 


Bodner  and  Partom  expressed  as  follows: 

*f\ 


(2.20) 
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where  DQ  is  the  limiting  value  of  plastic  strain  rate  in  shear,  2  is 
the  measure  of  material  hardening,  and  n  is  a  constant  controlling 
strain  rate  sensitivity. 

Being  strain  rate  sensitive  and  independent  of  a  specific  yield 
surface,  the  Bodner-Partom  flow  law  allows  plastic  deformation  during 
unloading.  This  is  a  significant  departure  from  classical  plasticity 
theory  where  all  unload  sequences  are  accomplished  elastically.  The 
strain  rate  sensitivity  exponent  (n)  influences  the  level  of  the  stress- 
strain  curves.  Decreasing  the  value  of  n  results  in  increased  strain 
rate  sensitivity  and  spreads  the  stress  strain  curves.  It  should  be 
pointed  out  that  the  strain  rate  sensitivity  exponent  (n)  is  a  material 
constant  and  independent  of  load  history  [4]. 

Thus,  the  Bodner-Partom  model,  henceforth  called  the  Bodner  rlfodel, 
accounts  for  both  isotropic  and  kinematic  hardening  and  is  capable  of 
representing  cyclic  load  effects. 

Bodner 's  Z  hardness  parameter  can  be  interpreted  as  a  macroscopic 
hardening  function  which  controls  resistance  to  plastic  flow.  Being 
deformation  history  dependent,  the  Z  parameter  is  assumed  to  be  a 
function  of  plastic  work: 

Z  -  Z(Wp)  (2.21) 

where  Wp  is  the  relative  amount  of  plastic  work  done  from  some  initial 
"•  state.  This  work  function  takes  the  form: 

Wp  ■  J  Sij'lPjdt  ■  2<b2j2>*  (2-22) 
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The  hardness  parameter  Z  now  becomes 

Z  •  Zx  -  (ZrZ0)exp  [-mWp]  (2.23) 

where  Z^  is  the  maximum  value  of  hardness,  ZQ  represents  the  initial 
value  of  hardness  form  which  plastic  work  is  measured,  and  m  is  the 
hardening  rate  exponent. 

Most  materials  exhibit  some  thermal  recovery  of  hardening  or  relax¬ 
ation  of  accumulated  plastic  work  at  high  temperatures.  In  order  to 
model  this  behavior,  the  form  of  plastic  work  must  be  redefined  as 

WP  ‘  |  sij^dt  + 1  fepry dt  (2-24> 

where 

2rec  '  -a(t7)  zl  <2'25> 

Z2  is  the  stable  non-work  hardened  value  of  Z  at  a  given  tempera¬ 
ture.  A  and  r  are  material  constants  picked  to  match  creep  test  data. 
The  thermal  hardening  recovery  term  (2rec)  of  Eq.  2.25  is  always 
negative  due  to  the  negative  sign  on  A,  since  Z  is  always  greater  than 
or  equal  to  Zg.  A  balanced  condition  exists  when  the  rate  of  work 
hardening  equals  the  rate  of  hardening  recovery.  At  high  strain  rates, 
Z  hardness  reaches  its  saturation  value  Z ^  quickly  and  a  steady  state 
condition  is  realized. 

Stouffer  experimentally  developed  the  nine  (9)  Bodner  Material 
parameters  for  IN-100  at  1350°F  [6]  and  are.  presented  in  Table  2.1. 
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TABLE  2.1 

BODNER  COEFFICIENTS  FOR  IN-lOO  AT  1350°F 


Material 

Parameter 

Description 

Value 

E 

Elastic  modulus 

26.3xl02KSI 

(18.133xl0^MPa) 

n 

Strain  rate  exponent 

0.7 

Do 

Limiting  value  of  strain  rate 

A 

10  sec 

zo 

Limiting  value  of  hardness 

915.0KSI  (6304  MPa) 

Z1 

Maximum  value  of  hardness 

1015. 0KSI  (6993  MPa) 

h 

Minimum  value  of  hardness 

600. OKS I  (4134  MPa) 

tn 

Hardening  rate  exponent 

2.57KSI'1  (.37273 
MPa-1) 

A 

Hardening  recovery  coefficient 

1.9x10" 3  ^c"1 

r 

Hardening  recovery  exponent 

2.66 

(1  KBAR  *  .100  MPa  *  14.504KSI) 

III.  METHODS  OF  ANALYSIS 


The  Computer  Program 

An  analytical  computer  program  of  known  accuracy  named  VISCO  was 
used  throughout  this  study.  VISCO  is  a  two-dimensional  finite  element 
code  employing  constant  strain  triangular  elements  for  both  plane  stress 
and  plane  strain  solutions. 

The  Bodner-Partom  viscoplastic  constitutive  equations  are  solved 
using  the  Gauss-Seidel  iterative  equation  solver  with  overrelaxation, 
eliminating  costly  stiffness  matrix  factorization  and  allowing  terms  of 
the  stiffness  matrix  to  be  changed  as  nodes  are  released  simulating 
crack  growth.  VISCO  employs  the  residual  force  technique  using  an 
elastic  stiffness  matrix  for  the  entire  analysis.  Plasticity  is  treat¬ 
ed  as  an  applied  load  used  in  conjunction  with  thermal  and  applied 
mechanical  loads  for  equilibrium.  During  each  time  step,  equilibrium 
tolerances  are  checked.  If  the  tolerances  are  exceeded,  the  time  step 
is  reduced  until  equilibrium  Is  obtained.  Time  integration  of  the 
Bodner  equations  is  accomplished  for  each  element  using  a  Euler  extrap¬ 
olation  scheme  [7]. 

To  accurately  model  cyclic  material  behavior  within  an  aircraft 
gas-turbine  engine,  it  is  necessary  to  understand  the  duty  cycle  of  the 
engine.  The  duty  cycle  can  be  characterized  by  frequent  load  cycles 
and  long  dwell  times.  In  addition,  few  rotating  jet  engines  components 
are  subjected  to  compressive  loading.  With  this  in  mind,  the  engine 
cyclic  load  spectrum  can  be  represented  as  a  saw-toothed  constant  am¬ 
plitude  load-time  pattern  with  non-zero  mean  stress  and  a  positive 
stress  ratio  as  shown  In  Fig.  3.1. 


TIME 


Fig  3.1  Typical  Cyclic  Load 

VISCO's  loading  function  was  modified  for  cyclic  response  (see 
Appendix  A).  Using  a  ramp  function  to  model  the  load  spectrum  in  Fig. 
3.1,  a  percentage  of  total  load  was  calculated  at  each  timestep. 
Element  stress,  strain,  and  Bodner  Z  hardness  were  recovered  at  each 
load  percentage.  The  load  rates  were  matched  to  the  half-cycle  period 
of  the  saw-toothed  load  function.  VISCO  normally  loads  the  specimen 
at  the  rate  of  20  percent  of  total  load  per  second.  In  order  to  match 
the  cyclic  load  spectra,  load  rates  were  matched  to  the  cyclic  half 
period.  For  example,  at  2.5Hz,  the  cyclic  half  period  is  .2  seconds. 
The  load  varies  from  a  minimum  to  peak  value  in  .2  seconds,  resulting 
In  a  500%  load  rate.  This  technique  graphically  depicts  incremental 
fatigue  effects  as  well  as  the  change  in  stress  field  at  a  crack  tip. 
In  addition,  modifications  to  VISCO  were  made  to  stop  load  cycling  at 
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a  predetermined  time  and  hold  the  load  at  a  constant  stress  intensity 
factor  (K). 

Computer  Program  Verification  and  Uniaxial  Studies 

To  verify  proper  program  operation  and  to  determine  IN-100  uniaxial 
non-linear  response  to  cyclic  loading,  a  two  triangular  finite  element 
model,  shown  in  Fig  3.2,  was  used  in  VISCO.  Plane  stress  solutions 
were  obtained  for  each  frequency  and  load  shown  in  Table  3.1. 


Fig  3.2  Uniaxial  Finite  Element  Model 


TABLE  3.1 


UNIAXIAL  TEST  LOADS  AND  FREQUENCIES 


Test  Frequency  (Hz' 

Stress ( ps 

2.50 

130,000 

0.167 

150,000 

0.03 

165,000 

0.03 

180,000 

0.03 

200,000 

0.03 

220,000 
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Peak  load  stress  levels  represent  typical  values  expected  at  the 
crack  tip  for  a  compact  tension  specimen  while  the  lowest  load/stress 
level  (130KSI)  represents  the  experimental  yield  stress  of  IN-100  at 


732°C. 

Each  of  the  eighteen  uniaxial  cases  was  allowed  to  run  for  fifteen 
complete  load  cycles.  Results  showed  that  the  computer  solutions  fol¬ 
lowed  the  input  load  cycle  at  each  time  step  and  that  calculated  stress 
values  were  accurate.  Plastic  deformation  calculations  compared  favor¬ 
ably  with  predicted  values. 

Computer  Computational  Time 

It  was  originally  intended  to  match  laboratory  fatigue  tests  done 
by  Larsen  and  Nicholas  with  computer  generated  fatigue  data  [14].  ;  ,s 

necessitated  at  least  8,000  load  cycles  before  goin*;  into  a  stained 
hold  period.  Preliminary  computer  runs  using  the  finite  element  mesh 
shown  in  Fig.  3.4  at  a  frequency  of  2.5Hz  required  2,000  seconds  central 
processor  time  per  run  on  the  CDC  6600  computer  for  5  complete  load 
cycles.  With  optimum  timestepping,  it  would  take  14  to  18  hours  to 
accomplish  8,000  load  cycles  of  central  processor  time.  Increasing  the 
size  of  the  finite  element  mesh  {reducing  the  number  of  elements)  doubled 
the  number  of  load  cycles  possible,  but  could  not  provide  critical  stress, 
strain  and  Z  hardness  values  near  the  crack  tip. 

With  these  computational  restrictions  in  mind,  it  was  decided  to 
model  the  cyclic  fatigue  effects  through  4  load  cycles  and  compare  re¬ 
sults  with  uniaxial  computations.  Changes  in  the  stress,  strain,  and 
Z  hardness  could  be  closely  monitored  with  each  cycle  and  changes 
observed. 


Two-Dimensional  Compact  Tension  Finite  Element  Modeling 

Two-dimensional,  fatigue  modeling  was  accomplished  using  a  standard 
compact  tension  specimen  geometry  as  shown  in  Fig.  3.3. 

Due  to  symmetry,  only  half  of  the  compact  tnesion  specimen  was 
modeled  using  constant  strain  triangular  elements.  The  finite  element 
mesh  (543  elements  and  327  nodes)  is  shown  in  Fig's.  3.4  and  3.5  and  is 
identical  to  mesh  3  used  by  Smail  for  compact  tension  specimens  of  the 
same  geometry  [16].  This  pattern  allows  unlimited  element  size  reduc¬ 
tions  and  insures  that  no  two  adjacent  elements  differ  in  size  by  more 
than  a  factor  of  2  [8].  Except  for  elements  near  the  loading  pin  holes, 
element  aspect  ratios  varied  from  0.5  to  1.0.  Elements  near  the  crack 
tip  had  an  area  of  4.8848xl0“®in^  (3. 1494x10" 5cm^ ) . 

Since  acceptable  accuracy  validation  was  performed  by  Smail  [15] 
using  the  same  computer  program  (VISCO)  and  finite  element  mesh,  no 
extensive  accuracy  tests  were  carried  out  herein.  Five  cyclic  load 
cases,  with  specified  stress  intensity  factors  shown  in  Fig.  3.6,  were 
studied.  Cyclic  stress  levels  at  the  crack  tip  were  compared  with  mono¬ 
tonic  computer  runs  made  by  Smail  at  the  same  load  and  stress  intensity 
factor  (see  Table  3.2). 
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STANDARD  COMPACT  TENSION  SPECIMEN 

Fig  3.3  Compact  Tension  Specimen  Geometry 


BASELINE  NO  DAMAGE 
K=35  KSI  t/in" 


R=O.I  FOR  ALL  CYCLING 


2.5  HZ 

K=35  KSI 


2.5  HZ 

K=45  KSI-y/lN 


.  167  HZ 

K=35KSI-yiN 


.03  HZ 

K=  35  KSI-y/lN 


Fig  3.6  Cyclic  Load  Cases  and  Stress  Intensity  Factors 
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IV.  RESULTS  AND  DISCUSSION 

At  2.5Hz  with  stress  levels  at  or  below  165KSI  (see  Fig's.  4. 1-4. 3), 
the  uniaxial  model  behaved  elasticity.  After  approximately  2.5  cycles, 
the  stress  strain  curve  loops  became  coincident.  A  quasi-stable  steady 
state  had  been  reached  in  as  little  as  2.5  load  cycles.  Stress  levels 
above  165KSI  caused  rapid  plastic  deformation,  especially  during  the 
first  load  cycle.  Subsequent  load  cycles  produced  constant  amounts  of 
inelastic  straining  as  shown  in  Fig's.  4. 4-4. 6.  No  stable  solution 
existed  for  these  stress  levels. 

As  load  cycle  frequency  decreased  no  stable-elastic  material 
behavior  was  noted  for  stress  levels  above  130KSI.  See  Fig's.  4.7-4.18. 
However,  for  all  cases  above,  the  greatest  plastic  deformation  always 
occurred  during  the  first  load  cycle.  Subsequent  cycling  produced 
constant  plastic  deformation. 

All  stress-strain  curves  displayed  linear  elastic  behavior  below 
130KSI  stress  levels.  Plastic  flow  appeared  to  start  near  the  material 
yield  stress  (130KSI  nominal  value)  in  all  cases.  Using  this  value  as 
a  reference  yield  stress,  specimens  cycled  at  .03Hz  spent  82  times 
longer  above  the  reference  value  than  the  specimen  at  2.5Hz.  Thus, 
it  was  expected  that  time  dependent  inelastic  deformation  was  dependent 
on  the  inverse  of  the  load  cycle  frequency. 

As  shown  in  Fig's.  4.19-4.21,  Z  hardness  is  very  sensitive  to  both 
frequency  and  load  level.  At  2.5Hz,  Z  hardness  increased  slowly  at 
stress  levels  below  165KSI,  while  rather  rapid  increases  were  noted  at 
lower  frequencies.  Increasing  the  stress  level  caused  rapid  rises  in  Z 
hardness.  Z  hardness  saturation  occurred  after  3  cycles  in  all  cases  for 
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stress  levels  of  180KSI  or  greater,  with  the  exception  of  the  2.5Hz 
specimen  at  180KSI - 

If  one  compares  the  cyclic  stress-strain  loops  of  Fig's.  4.1-4.18 
to  the  Z  hardness  curves.  Fig's.  4.19-4.21  show  constant  magnitude 
plastic  strain  per  cycle  occurs  only  when  Z  hardness  saturates. 

In  addition,  Z  hardness  does  not  oscillate  with  the  load  cycle. 

The  plastic  deformation  occurring  during  each  load  cycle  causes  Z  to 
increase  until  reaching  saturation.  After  Z  reaches  saturation,  the 
additional  plastic  deformation  accumulated  per  cycle  is  balanced  by 
the  relaxation  due  to  temperature.  Thus,  the  material  has  reached  its 
fully  hardened  state  and  does  not  vary  with  additional  plastic  deforma¬ 
tion.  The  foregoing  results  graphically  depict  the  thermal  hardening 
recovery  and  secondary  creep  characteristics  of  the  Bodner  model.  By 
setting  the  time  derivative  of  equation  2.24  equal  to  zero,  the  rate  of 
work  hardening  equals  the  rate  of  thermal  hardness  recovery.  For  a  given 
time  increment  after  L  reaches  saturation,  a  balanced  condition  is 
present  which  represents  secondary  creep.  Each  additional  load  cycle 
causes  constant  magnitude  plastic  deformation  increases  per  cycle, 
shifting  the  stress-strain  curves  to  the  right  (Fig  4.1  -  4.18). 

Comparing  stress  levels  and  time  values  in  Fig's  4.19  -  4.21  reveals 
that  7  saturates  at  approximately  the  same  time  and  is  independent  of 
frequency. 
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It  should  be  noted  in  Fig.  4.22,  which  is  typical  of  all  uniaxial 
cases,  that  plastic  strain  increases  during  that  portion  of  the  loac4 
cycle  where  the  stress  is  above  the  reference  yield  stress.  These 
strain  rise  points  occur  at  approximately  80  percent  of  the  applied 
load.  During  the  unload  portion  of  the  cycle  while  above  the  yield 
stress,  plastic  deformation  continues.  Between  cycles,  plastic  strain 
remains  constant  until  the  subsequent  "yield  point"  is  reached.  The 
magnitude  of  the  rise  in  plastic  strain  during  each  load  cycle  de¬ 
creases  until  Z  hardness  reaches  its  saturation  value  (maximum). 
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Fig  4.22  Plastic  Strain  vs  Time 


Compact  Tension  Specimen  Tests 

Results  and  comparisons  made  in  this  section  are  presented  for 
each  frequency  and  stress  intensity  factor  shown  in  Fig.  3.6.  Sub¬ 
sequent  figures  are  in  five  groups  as  shown  below: 

a)  Cyclic  stress  strain  curves  at  the  crack  tip  for 
the  first  3  load  cycles 

b)  Profiles  of  y  stress,  y  strain  and  2  hardness  as  a 
function  of  distance  ahead  of  the  crack  tip  for  the 
the  first  3  cycles 

c)  Changes  in  plastic  zone  size  with  the  number  of 
load  cycles 

d)  Vertical  (y)  displacement  behind  the  crack  tip  versus 
distance  behind  crack  tip  during  cyclic  unload 

e)  Cyclic  crack  mouth  displacement  versus  load  curves 

As  shown  in  Fig's.  4.23,  4.24,  4.25,  and  4.26,  (effective  stress 
versus  total  strain  curves);  rapid  increases  in  plastic  deformation  at 
the  crack  tip  occur  during  the  first  load  cycle.  No  significant  changes 
occur  afterward.  Constant  magnitude  plastic  deformation  occurs  with 
each  successive  cycle.  If  one  compares  these  results  with  the  Z  hard¬ 
ness  profiles,  Fig's.  4.29,  4.32,  and  4.38,  it  can  be  seen  that  the  Z 
hardness  near  the  crack  tip  reaches  saturation  in  the  first  half  load 
^  cycle.  The  material  at  the  crack  tip  is  fully  hardened;  further 

cycling  produces  no  further  hardening.  Referring  to  Fig's.  4.29-4.38, 
no  significant  changes  in  y  strain,  Z  hardness,  or  y  stress  profiles 
('■%  ahead  of  the  crack  tip  occur  after  the  first  cycle.  This  behavior 
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suggests  a  stable  condition  exists  at  the  crack  tip  after  the  first  cy¬ 
cle  and  can  be  compared  directly  with  uniaxial  results  shown  in  Fig's. 
4.1-4.18.  This  suggests  that  the  material  ahead  of  the  crack  tip  is 
essentially  in  a  stress  controlled  boundary  condition,  and  accumulates 
plastic  strain  per  cycle  just  as  it  does  in  the  uniaxial  cases.  No  net 
accumulation  of  strain  per  cycle  would  be  seen  if  the  material  ahead  of 
the  crack  tip  was  strain  controlled.  In  addition,  this  type  of  behavior 
ahead  of  the  crack  tip  infers  that  if  nodes  are  released  in  the  finite 
element  model  to  simulate  crack  growth,  and  the  number  of  load  cycles 
between  node  releases  is  large;  changes  ahead  of  the  crack  tip  will 
occur  the  instant  the  node  is  released.  One  cycle  after  node  release, 
the  stress-strain  field  and  Z  hardness  remain  essentially  constant  until 
the  next  node  is  released.  This  suggests  that  the  behavior  ahead  of  the 
crack  associated  with  crack  growth  can  be  characterized  by  a  short  tran¬ 
sient  period  after  node  release  followed  by  stable  cyclic  behavior  until 
the  next  node  is  released. 

If  the  stress  profiles  ahead  of  the  crack  tip  are  known,  one  can 
refer  to  the  uniaxial  results  and  come  very  close  to  predicting  the  be¬ 
havior  in  the  compact  tension  finite  element  model.  However,  the  stress 
ahead  of  the  crack  tip  using  the  Linear  Elastic  Fracture  Mechanics  (LEFM) 
plane  stress  solution  of  Eq.  4.1,  [1,11] 

+  sin  cos  ■yj  (4.1) 
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with  theta  (e)  of  zero  degrees,  stress  intensity  factors  (K^)  of  35  and 
45KSIvTn  respectively  differed  significantly  from  the  finite  element 
solution.  Results  depicted  in  Fig's.  4.27,  4.30,  4.33,  and  4.36  show 
that  the  LEFM  solution  overestimates  the  stress  field  near  the  crack 
tip  in  most  cases  until  reaching  0.02  inches  ahead  of  the  crack  tip. 
Beyond  0.02  inches,  LEFM  underestimates  stresses  by  10%.  Very  near  the 
crack  tip  (under  0.006  inches),  the  LEFM  solution  overestimates  the 
stress  by  approximately  15  percent. 
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Fig  4.25  Effective  Stress  vs  Total  Strain 
at  Crack  Tip  .167Hz  K=35KSI/Tn 
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Results  of  a  typical  overload  condition  are  displayed  in  Fig's.  4.39- 
4.41.  The  load  was  cycled  at  a  stress  intensity  factor  of  45K$l/fn  for 
two  cycles.  Load  cycling  was  stopped  when  the  stress  intensity  factor 
rose  to  35KSI/Tn  during  the  next  half  cycle.  Stress,  strain  and  Z  hard¬ 
ness  values  shown  represent  the  effect  of  overload.  It  can  be  seen  that 
overload  increases  plastic  deformation  near  the  crack  tip,  retarding 
future  crack  growth.  This  verifies  laboratory  overload  fatigue  crack 
growth  studies  conducted  by  Larsen  and  Nicholas  [14]. 

Numerous  studies  have  shown  crack  closure  in  center  cracked  panels 
during  cycles  unload  [15].  To  check  for  crack  closure,  vertical  (y) 
displacements  behind  the  crack  tip  were  continuously  monitored  during 
the  unload  cycle.  As  shown  in  Fig's.  4.42-4.46  no  negative  displacements 
were  observed,  indicating  no  crack  closure.  Lack  of  closure  may  be 
attributed  to  both  specimen  geometry  and  to  the  positive  load  levels 
(stress  ratio  of  0.1)  throughout  the  unload  half  cycle.  The  compact 
tension  specimen  geometry  can  be  thought  of  as  a  cantilever  beam  with  a 
positive  end  load  with  the  fixed  support  at  the  crack  tip.  With  .this 
in  mind  It  can  be  seen  that  one  would  not  expect  closure  (negative) 
displacements  at  the  crack  tip. 
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Fig  4.40  Y  Strain  vs  Distance  Ahead  of  Crack  Tip 
Considering  Overstress 
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Fig  4.42  Y  Displacement  vs  Distance  Behind  Crack  Tip 

During  Unload 
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Throughout  this  study,  it  was  speculated  that  at  stress  levels 
below  165KSI  near  the  crack  tip,  the  compact  tension  specimen  behaved 
elastically  and  that  plastic  deformation  produced  negligible  effects. 

To  verify  this  premise,  an  averaging  process  was  undertaken.  The 
load  was  cycled  through  4.5  cycles.  After  4.5  cycles,  the  finite 
elements  which  had  plastic  strain  values  at  or  above  10  in/in  were 
considered  plastic.  By  averaging  stresses  and  total  strains  for  all 
plastic  elements,  the  stress-strain  curves  (Fig's.  4.46-4.47)  were 
plotted.  At  2.5H z  and  Ka35KSI  in  the  compact  tension  specimen  reached 
a  "quasi-steady  state"  after  2.5  cycles.  Additional  cycling  produced 
negligible  shifts  in  the  stress-strain  curve.  Peak  stresses  averaged 
165KSI  at  the  crack  tip  over  the  plastic  zone.  At  0.03Hz,  no  such 
"steady  state"  condition  was  evident.  Similar  observations  were  made 
with  uniaxial  results  at  the  same  stress  levels. 

To  observe  overall  effects  of  load  cycling,  extended  hold,  and 
node  release,  the  load  was  cycled  at  2.5Hz  at  a  stress  intensity  factor 
of  35KSI  in  for  2.5  cycles.  After  cycling,  the  load  was  held  at  a  con¬ 
stant  stress  intensity  factor,  for  fifteen  minutes  followed  by  node 
release.  Results  shown  in  Fig's.  4.48-4.49  depict  plastic  zone  growth 
with  the  load  cycling.  Comparing  Fig’s.  4.49  with  4.50  shows  there  is 
little  change  in  plastic  zone  size  during  the  sustained  hold  period. 


AVERAGE  TOTAL 


The  stress  field  ahead  of  the  crack  tip  changes  little  during  the 
sustained  hold  and  during  node  release  (see  Fig.  4.52).  However,  strain 
ing  continues  to  increase  slightly  (Fig.  4.53).  It  should  be  noted  that 
after  node  release,  the  Y  strain  (e^)  is  coincident  with  the  strain 
curve  prior  to  node  release.  In  addition,  it  must  be  realized  that  the 
stress-strain  field  has  translated  one  finite  element  to  the  new  crack 
tip.  The  plastic  zone  ahead  of  the  crack  tip  is  approximately  the  same 
size  as  it  was  prior  to  node  release.  (Fig's.  4.50-4.51)  As  expected, 
the  Z  hardness  parameter  decreased  during  hold  period  (see  Fig.  4.54). 
During  the  sustained  hold,  the  rate  of  plastic  deformation  (strain  rate) 
decreased  due  to  the  absence  of  load  cycling,  allowing  the  rate  of 
thermal  hardness  recovery  (2  ,  equation  2.24)  to  become  more  dominant. 

ICv 

The  decrease  in  Z  hardness  (strain  hardening)  indicates  a  gradual 
return  toward  an  unhardened  state. 

Changes  in  plastic  zone  size  with  frequency  are  shown  in  Fig's. 
4.48-4.49  and  4.55-4.60.  It  should  be  noted  that  as  the  frequency  de¬ 
creased  the  size  of  the  plastic  zone  increased.  This  corresponds 
directly  to  uniaxial  test  results  shown  in  this  section. 

Laboratory  fatigue  testing  relies  heavily  on  compliance  measure¬ 
ments  to  acquire  data.  Little  information  is  available  as  to  compact 
tension  specimen  far  field  behavior.  The  far  field  behavior  dictates 
the  linearity  or  non-linearity  of  clip  gauge  measurements  taken  at  the 
crack  mouth.  To  this  end,  computer  generated 'far  field  displacements 
(Y  displacement)  at  the  crack  mouth  (node  7)  on  the  edge  of  the  speci¬ 
men  were  plotted  during  the  load  cycle  for  three  frequencies  and  two 
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stress  intensity  factors.  The  results  show  (Fig's.  4.61-4.64)  that  the 
compact  tension  specimen  behaves  elastically  at  points  on  the  specimen 
boundary  and  does  not  see  the  localized  inelastic  behavior  near  the 
crack  tip. 

To  check  the  finite  element  compliance  data,  elastic  modulus  (E) 
calculations  were  made  using  compliance  data  from  Fig's.  4.61-4.64,  and 
the  load  and  crack  length  data  from  Table  3.2  with  the  experimentally 
derived  equation  (5.1)  for  elastic  compliance  calculations  for  a  compact 
tension  specimen  [16]: 


27.2 


E 


(5.1) 


where  B  is  the  specimen  thickness,  C  is  the  compliance  and  E  is  the 
elastic  modulus  for  the  material.  Using  the  specimen  thickness  of 
.2154  inches  and  compliance  of  4.4842xl0’6  in/lb,  the  calculated  modu 
lus  was  28.1xl06  PSi.  By  comparing  the  calculated  elastic  modulus  to 


the  experimental  elastic  modulus  (Table  2.1),  it  can  be  seen  that  the 
calculated  elastic  modulus  is  7  percent  greater  than  experimentally 
determined  elastic  modulus. 


Fig  4.61  Crack  Mouth  Displacement  vs  Percent  Load 
2.5Hz  K=45KSI/Tn 
Normalized  to  K=35KSI/TrT 
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Fig  4.64  Crack  Mouth  Displacement  vs  Percent  Load 


V.  CONCLUSIONS 


The  foregoing  computations  are  based  on  the  Bodner  constitutive 
equations  which  model  the  real  behavior  on  IN-100  at  elevated  temper¬ 
atures  in  which  experimental  tests  have  been  made.  Computations  show 
that  this  material  has  a  significant  amount  of  time  dependent  inelastic 
behavior  as  evidenced  by  computational  results  for  both  uniaxial  and 
compact  tension  specimens. 

1.  There  is  close  correlation  between  uniaxial  cyclic  behavior 
and  the  material  behavior  near  the  crack  tip. 

2.  At  frequencies  above  .167Hz  and  at  stress  levels  below 
165KSI ,  inelastic  behavior  is  negligible  with  respect  to 
total  material  behavior  for  both  uniaxial  and  compact 
tension  tests. 

3.  At  high  stress  levels, (independent  of  frequency)  large 
amounts  of  plastic  deformation  occur  in  the  first  load 
cycle.  Subsequent  load  cycles  produce  constant  amounts  of 
plastic  deformation  per  cycle. 

4.  The  stress  field  and  plastic  zone  size  ahead  of  the  crack 
tip  remains  relatively  constant  after  one.tothree  load 
cycles  in  the  compact  tension  geometry. 

5.  Linear  Elastic  Fracture  Mechanics  overestimates  the  stress 
field  directly  ahead  of  the  crack  tip  and  underestimates 
it  at  larger  distances. 
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6.  The  stress  field  ahead  of  the  crack  tip  becomes  more 
uniform  during  sustained  hold  periods. 

7.  The  size  of  the  plastic  zone  remains  relatively  constant 
during  sustained  hold  periods. 

8.  After  node  release  from  a  sustained  hold  period,  the 
size  of  the  plastic  zone  ahead  of  a  crack  tip  remains 
essentially  constant  and  only  shifts  forward  to  the  new 
crack  tip. 

9.  The  size  of  the  plastic  zone  after  cycling  is  frequency 
and  load  dependent.  For  a  given  load,  the  plastic  zone 
increases  in  size  as  the  frequency  decreases. 

10.  The  cyclic  behavior  near  the  crack  tip  for  compact  tension 
geometry  is  stress  controlled. 

11.  The  far  field  overall  behavior  of  the  compact  tension 
specimen  is  elastic  and  independent  of  frequency  and 
load  level. 
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Appendix  A 

The  computer  program  VISCO  was  modified  for  cyclic  loading  by 
adding  a  subroutine  named  CYCLIC  and  appropriate  call  statements  in 
the  Load  Subroutine.  CYCLIC  uses  a  constant  sloped  ramp  to  simulate 
the  input.  The  slope  of  the  ramp  is  computed  using  the  half-cycle 
period  of  the  input  frequency.  A  local  time  variable  is  used  to  lo¬ 
cate  position  on  the  input  cycle.  At  the  peak  load  the  slope  is  re¬ 
versed  and  the  local  time  frame  shifted  to  the  input  peak.  During 
each  tlmestep,  (passed  from  the  main  program)  subroutine  CYCLIC  cal¬ 
culates  the  percentage  of  total  load.  The  calculated  load  percentage 
is  then  passed  back  to  the  Load  Subroutine.  Nodal  loads  were  computed 
and  used  to  calculate  appropriate  displacements,  and  stresses  and 
strains.  The  process  is  repeated  for  each  timestep. 

Variables  read  Into  the  program  in  statements  labeled  one  and  two 
control  load  cycling  and  stop  times. 

The  variable  STOPCY  is  the  time  variable  used  to  stop  load  cycling. 
NC  dictates  whether  subroutine  CYCLIC  Is  called  from  the  Load  Subroutine 
(statement  labeled  3).  If  NC  equals  one,  cyclic  loads  are  applied.  Any 
other  number  causes  a  monotonic  load  to  be  applied.  Variables  PMAX,  PO, 


and  PERIOD  control  the  Input  load  cycle.  PMAX  is  the  peak  load  percent 
age  (100%),  while  PO  is  the  stress  ratio  for  the  load  cycle.  The 
variable  PERIOD  Is  the  cycle  period  of  the  Input  frequency  in  seconds. 
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